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ABSTRACT 

We examine the two-point correlation function (TPCF) of local maxima in tempera- 
ture fluctuations at the last scattering surface when this stochastic field is modified by 
the additional fluctuations produced by straight cosmic strings via the Kaiser-Stebbins 
effect. We demonstrate that one can detect the imprint of cosmic strings with tension 
Gfi > 1.2 x 10 -8 on noiseless 1' resolution CMB maps at 95% confidence interval. 
Including the effects of foregrounds and anticipated systematic errors increases the 
lower bound to Gfi > 9.0 x 10 -8 at 2a confidence level. Smearing by beams of order 
4' degrades the bound further to Gfi > 1.6 x 10 -7 . Our results indicate that 2-point 
statistics are more powerful than 1-point statistics (e.g. number counts) for identifying 
the non-Gaussianity in the CMB due to straight cosmic strings. 

Key words: Cosmic Microwave Background, two-point correlation function, cosmic 
strings, peaks, non-Gaussianity 



1 INTRODUCTION 

The origins of the seeds of present large scale structures 
in the universe are still debated. It is believed that they 
are mainly primordial and produced some time after Big- 
Bang. In this framework, there are two approaches: 1) the 
freezing- in of quantum fluctuati ons of a scalar field during 
the s o -called inflationary epoch ( Guthlll98ll ; I Liddle fe Lvthl 
Il993l : I Liddld [l999l : I Steinhardl 19951): and/or 2) top olog- 
ical defects as sources (|Kibblel Il976l : I Kibble) Il980h . In- 
deed topological defects can be formed during phase tran- 
sitions between different vacuum states in an expand- 
ing universe, and cosmic strings are predicted by quan- 
tum field theory in cosmology | Bevis et al. 120081 pOlOj : 



DePiesll2009l:lHindmarsIi fe Kibbldl995l:ISakellariadoij| 2006 
Vachaspati fe Vilenkin 1984; Shellard 1987; Khlopov 1999 



Vilenkin fe Shellardl l200ol: iKibbld 1 19761 : IZeldovichl [l98L 
Vilenkinll 1 98 ll . 1 1 9851 : [Allenll 1 9971 ) . Both the inflationary and 
topological defects scenarios predict the same features for 
cosmic microwave background (CMB) power spectrum on 
large scales. But at the intermediate and small scales, due 
to differences in the super-horizon scale behavior of per- 
turbations in these theories, the predictions are completely 
different. 

The inflationary ACDM paradigm is consistent with 
today's high precision observations of the CMB. Neverthe- 
less, from both theoretical and observational points of view 
there are many motivations for other sources of anisotropics. 



in hybrid inflation models, the brane- world paradigm 
superstring theory, producti on of topologic al de- 



E.g., 
and 

fects are crucial and inevitable (I Copeland et al.l 11994 



Sarangi fe Tve 2002; Tye 2008; Kibble 2004; Copeland et al. 



| 2004l: iMaiumdar fe Davisl 120021: IPogosian et al.l 120031 : 
iDvali fe Vilenkinll2004l : ISakellariadoJl997ll . For more recent 
observational results see^evis et al.l (|2Q08h ; I Dvorkin et al.l 
(|201lh : lRingeval fe Bouchetl <|2012h . 



A cosmic string (CS) network which consists of infinite 
strings, loops and junctions of strings can generate gravita- 
tional waves as the universe evolves. Astrophysical evidence 
of CS depends on 1) the inter-commuting probability and 2) 
the dimensionless string tension Gfi/c 2 = A 2 /Mp lanck , where 
G is Newton's constant, \i is the mass per unit length of the 
cosmic string and A is the energy scale of the string creatio n 
epoch (| Bevis et al.ll2008l l2010l : IVilenkin fe Shellardl l2000h . 
We set c = 1 throughout this paper. Determining bounds 
on the value of \i directly means limiting the basis of fun- 
damental theories for CS production. In addition, observing 
cosmic strings not only is a kind of observational evidence 
for such theories, but also provides an opportunity to rule 
out or confirm theoretical models of particle physics. 

There are many constraints on the upper as well as 
lower values of cosmic string parameters from theoreti- 
cal and observational perspectives. Pulsar timing and pho- 
tometry based on gravitational microlensing and gravita- 
tional waves require 10" 15 < Gfi < 10" 8 (Jenet et al.l 
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20061: Pshirkov & Tuntsov 2010: Tuntsov 


& Pshirkovll2010l: 


Damour & Vilenkin 2005; Battve & Moss 


201(4 Oknvanskiil 


1999:lGasilov et aljll985: Blinnikov et all 


1982). The COS- 



MPS survey requires Gfi < 3 x 10 7 ([Christiansen et aD 
2010 ). The 21cm signature of CS has been investigated in 
(|Branderberger et al.1 lioiol ). The LIGO and VIRGO col- 
lab or oXlon^hm^deteY^iie^l x 10 -9 < G\i < 1.5 x 
10" 7 (jLIGO and VIRGOl [jpoj ). In a very recent paper, 
the stochastic gravitational waves from the European pul- 
sar timing array p lace constrains G\i < 5.3 x 10 -7 
(jSanidas et alJl2012h . 

Another strong constraint on Gfi comes from the in- 
vestigation of temperature fluctua tions at the last scat- 
tering surface (|Fraisse et al.l 120081 ). The accumulation of 
anisotropics induced by CS on the fluctuations at last 
scattering surface can be divided into two categories: 1) 
anisotropics related to pre-recombin ation processes and cre - 
ated by the Kaiser- Stebbins effect ([Kaiser fe Stebbindfl984h 
and 2) the decay of string loops, which results in a stochas- 
tic background of gravi tational waves ( Frai sse" et al.l l2008: 
iKaiser &; Stebbins1ll984l ). CMB analyses bo und the cosmic 
strin g tension to be Gi l < 6.4 x 10~ 7 ([Battve Moss 



2010 



2005 



2005, 



Fraisse et aD 1200a IKaiser fe Steb bins 198 3: I Fraisse 
Bevis et al 1 120041: iBevis et al.1 ^PoT T Wyman et al 



2006: |PerivolaroDoulodll993aj ). The effects of CS 



the skewness of the one-point probability distribution of 
CMB tem peratures ([ Yamauchi etajj|2010aj), the TT power 
spectrum (Yamauchi et al. 2010b|), the B-mode polarization 
(|Ma Brownll2010l ) have also been considered. 

Various detection methods have been explored, in- 
cluding: Wav elet domain Bayesian denoising: G\i > 
6.3 x 10~ 10 (| Hammond et al J l2008 | ). The Canny algo- 
rithm: Gu > 5.5 x 10~ 8 (bangs fc Brandenbergertl2010al: 



Stewart &: Brandenbergerl 20091 : Danos &; Brandenbergerl 
201 Obi ). Level crossing analysis: G\i > 4 x 10~ 9 and G\i > 
5.8 x 10 -9 wi thout and in the presence of instrumental noise, 
respectively (|Movahed fe Khosravil liml ). 

More recent observational constraints via WMAP 
and the South Pole Telescope yield Gji < 1.7 x 
10~ 7 at 95% confidence interval (I Dvorkin et al.1 l201ll ). 
Other computations expressed G\i < 0.7 x 10 -6 with 
fe=io < 0.11 (the fractional contribution of cosmic 
string on th e temp erature power spectrum at £ = 10) 
(| Bevis et al.l 120081 ). The constrai nts on CS fro m fu- 
ture CMB polarization come from ([Foreman et al.l 1 20 111 ) . 
The non -Gaussianinty imposed by CS has been consid- 
ered bv (I Starck. Aghanim and Fornl 120041: iRingevaJ I20T0I : 
iHobson et al.l 1 19991 : iBarreiro and Hobsonl l200ll ). By using 
the WMAP 7-year data release and the SPT £ < 3000 
power spectrum an upper limit on cosmic string tension has 
been given as Gji < 1.7 x 10 -7 at 95% confidence interval 
(I Dvorkin et all Eoilh. In the same paper, the aspect of po- 
larization power spectrum to put robust constrains on the 
properties of CS have been discussed. 

In the current study we concentrate on the discontinu- 
ities and fluctuations in the CM B map arising by CS from 
the Kaiser-Stebbins (KS) effect (jKaiser fc Stebbind fl984l ) . 
This phenomenon can produce observational consequences 
on the anisotropics in the CMB with high degree of relia- 
bility in the high resolution map. In what follows we study 
the two-point correlation function of lo cal maxima or min- 
ima in the observed temperature maps ([Bond &; Efstathioul 



Il987l : iHeavens &; Shethll 19991 ) to see if this is a useful probe 
of the extra roughness in the temperature distribution in- 
duced by strings. Previous work has shown that although 
the bispectrum of all pixels in a map at 5.5' resolution is not 
sensitive to a CS component, the two-point correlation func- 
tion of local maxima is, especially on scales of order 10 — 20 
arc- minutes ([Heavens fe GuDtall200ll ). The main goal of the 
present work is to quantify the limits of G\i which such a 
measurement can place. 

Section 2 describes how we generate mock maps of pri- 
mordial Gaussian CMB, and how we incorporate the effects 
of straight cosmic strings. In essenc e this is a straightf o rward 
comb ination of the algorit h ms in IHeavens Shethl (| 19991 ) 
with iMovahed &; Khosravil (| 20 111 ) . We then identify peaks 
in these maps and study if the one- and two-point statistics 
of these peaks agree with the Gaussian prediction, quanti- 
fying the statistical significance of the differences. A final 
section summarizes why we conclude that 2-point statistics 
are much more efficient for identifying the presence of CS 
compared to 1-point statistics. 



2 SIMULATION AND ANALYSIS OF MOCK 
CMB MAPS 

This section describes how we simulate maps of th e last 
scattering surface (Tstewart Bran denberger 2009; 



Danos &: Brandenberge 
1993b; Moessner et al. 



2010 a|b|: | Per i volar op oulod 
~~ |Per i volar op oulod 1993al : 
Movahed &; Khosravil 1201 ll ) . At first, our code creates 
pure Gaussian fluctuations corresponding to the standard 
inflationary model with ACDM com p onents in a flat 
universe following iBond Efstathioul (| 19871 ). However, 
our program can be easily modified to other cosmo- 
logical models for this purpose. Secondly, anisotropics 
produced by straight cosmic strings by the Kaiser- 
Stebbins effect, from the last scatterin g surface up to 
the present, are simulated following IPeriyolaropoulosI 



(Il993al) and I Stewart fe Brandenbergerl d20Q9h: 
Ipanos &; Brandenbergerl (|2010al bl): I Movahed &; Khosravil 
( 2013 ). Our method of simulation differs from that used 
to pr oduce the CS maps analyzed by (Hea vens &; Guptal 
l200ll ). where the CS contribution to the ene rgy-momentum 
tenso r was modeled using Fourier methods ([Pogosian et al.l 
Il999h . 

For reasons di scusse d in detail in 

Stewart & Brandenberger (2009); iDanos h Branden bergerl 
(|2010al lbl): IMovahed fe Khosravil (|201ll ). to simulate the 
Kaiser-Stebbins lensing due to moving cosmic strings we 
work in real space, where straight strings produce random 
jumps on the background radiation field. The scaling be- 
havior of straight strings means that the number of strings 
crossing a given Hubble v olume is fixed to M s t r in g = 10 
([Bennett k Bouchetl [l988). The cosmic strings possess 
relativist ic velocities; consequently, after 2£h (£h is the 
Hubble time) an entirely new network of cosmic strings 
provides new kicks to the CMB photons. The two signals are 
superposed, then smeared by our model of the instrumental 
beam, after which we add instrumental noise. Finally, we 
identify peaks in the simulated maps and measure the 
peak-peak correlation, showing results after averaging over 
a large ensemble of realizations. 
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Figure 1. Contribution of various components to the total angu- 
lar power spectrum of temperature fluctuations. Solid line shows 
the power-spectrum for a WMAP-7 ACDM model (from CAMB). 
Symbols show a simulated Gaussian map with resolution R = 1' 
and size 10°. Long-dashed line shows the contribution from cos- 
mic strings having Gfi = 8 x 10 -8 ; dot-dashed line shows a power- 
_n qn+ - 05 

law t(t + 1)C £ ~ £ u - yu -o.o5. Dashed-dot-dot curve shows the 
window associated with a beam having FWHM= 20'. Clearly, 
the CS component is most easily detected at t ~ 6000. 



2.1 Mock Gaussian CMB map 

In what follows, the size and resolution of simulated map 
are G and i?, respectively. Thus, a G = 10° map at R — 1' 
requires 600 x 600 pixels. The rms instrumental noise, a no ise, 
and the full width half maximum of detector, FWHM, are 
used to take into account additional effects on the simulated 
maps. We set a no ise — which is appropriate for th e 

South Pole Telescope (jRuhl et al.ll2004l : iKeisler et al.ll201lh . 
We also use FWHM = 4' to illustrate our results. 

For making Gaussian maps, all that is required is the 
initial power spectrum C g. For this, we use the CAMB soft- 
ware ([Lewis et al.ll2000h with parameters appropriate for a 
ACDM model that is consistent with the WMAP-7, Su- 
pernova type la (SNIa) and the Sloan Digital Sky Survey 
(SDSS) datasets. We u se this Cg to generate a 2D Gaussian 
random field following iBond &; Efstathioul (|l987h . Since we 
are interested in relatively small angula r scales, we work 
i n the flat sky approximation, following Hea vens &; Sheth 
<|l999h . 

To add the effec ts of strings we follow 

(I Stewart &; Brandenberger 20091; iDanos &; Brandenbergerl 



l2010al lbl; iMovahed fe Khosravil 1201 lh . This means that we 
ignore the contribution of CS loops, since their size is 
smaller than our map resolution (1 arcmin). In contrast the 
characteristic length scale of straight CSs is the Horizon 
scale. 



2.2 Combination of simulated components 

When combining the Gaussian (G) and string (S) 
components, we are careful to ensure that, at £ < 
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Figure 2. Various components of simulated maps, for a WMAP- 
7 ACDM model, with G\x = 8 x 10 -7 strings. The map size is 
5° x5° at resolution R = 1' and smeared by a beam of FWHM=4'. 
Circles in some panels show peaks above # = 0.5cro- 



200 (IDanos fe Brandenbergerl lioiOal : IMovahed fc Khosravil 
l201ll ). the total power is the same as that observed. In prac- 
tice, this means that in each pixel (z,j) we set the fluctua- 
tion, T = (T(iJ) - (T»/(T), to be 



^(g+s) (h j) = uT {G) (i, j) + (i, j) 



(I) 



where uj < 1 is chosen so that the amplitude of the power 
spectrum 



cf +s) = ^c\ G) + c { t S) 



(2) 
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Figure 3. Abundance of density of peaks as a function of peak 
threshold level for pure Gaussian maps and upon adding a cos- 
mic string component. Upper left: G\x = 2 x 10 -8 . Upper right: 
GfJL = lx 1CT 7 . Middle left: Gfi = 8 x 1CT 7 . Results for a Gaus- 
sian map (Gaussian- (GS)) having the same power spectrum as 
the Gaussian+String map has also been indicated in this panel. 
Middle right: for Gfi = 8 x 10 -7 with beam effect. The lower 
panels show the difference between theoretical Gaussian field pre- 
diction for the number density of peak and what measured in the 
simulated maps. 

is close to that observed at i < 200. Since C\ S ^ depends on 
£, determination of the appropriate value of uj is done by a 
likelihood analysis. 

The beam smearing is modeled by: 



C (G+S)B = c (G+s) e -r 2 €(€+i) 



J t — c (3) 

with r FWHM/ V8ln2 ([Bond fe Efstathioul Il987l : 

iHeavens fe Shethlll999h . 

Finally, we add a model for the noise: 



(4) 



where the final noise term is white, i.e., it has 
(J r ( A r)(ri)J r ( i v)(r2)) ~ foirac(ri -r 2 ), with the noise in pixel 
being a zero-mean Gaussian number with rmS (Jnoise- 
Figure [2] illustrates various components and steps in our map 
making process. 



2.3 Peak counts in mock maps 

We have checked that the number density of peaks we iden- 
tify in the Gaussian maps agrees with that expected from 
theory. When the peak height is expressed in units of the 



rms temperature, & = AT/T, this predi ction depends only 
on th e shape of the power spectrum Gt ([Bond &; Efstathioul 
Il987h . Since CS modify Ci at high £, it is interesting to ask 
if the peak counts predicted by C[ G+S ^ provide a good de- 
scription of the peak abundances in the G + S maps, even 
though the maps themselves are not Gaussian. If not, then 
peak counts alone allow one to distinguish between a purely 
Gaussian model and one with an additional component. 

Figure [3] shows that, for G\i ^ 10 -7 (top panels), the 
peak counts in G and G+S are almost indistinguishable. For 
larger values of G/i, the peak counts are noticably different 
from one another (middle left panel), with the distribution 
being shifted to smaller mean values when CS are present. 
However, the measure ments are each well de scribed by 
Gaussian peaks theory ([Bond &; Efstathioulll987f ) with their 
respective power spectra (C[ G ^ for the circles, and C[ G+S ^ 
for the triangles), even though the G+S maps themselves are 
not Gaussian. Thus, given only the observed Ci and the peak 
counts, it will not be possible to determine if there is a CS 
component in the maps. Beam smearing makes the counts 
indistinguishable upto even larger values G\i ^ 8 x 10 ~ 7 
(middle right panel of Figure [3}. In addition, to make more 
obvious, we have added An($) = n CO m.{^) ~ nthe.{$) in the 
lower panel of Figure [3] Where "com" refers to numerical 
result and " th e." corresponds to the oretical prediction. 

Recently, IPogosvan et al.l ([201 if ) derived expressions for 
the number density of extrema in weakly non- Gaussian 2- 
Dimensional fields. They showed that various non-Gaussian 
models could be distinguished by means of n(#). Our anal- 
ysis demonstrates that, at least for the non-Gaussianity due 
to straight CSs, this does not work. 



2.4 Two-point statistics 

Although we have demonstrated that peak counts in our 
G + S maps are consistent with those in a Gaussian field 
having the same Ce, direct inspection of the maps themselves 
(top panel of Figure 2} shows that they have quite different 
morphologies. The CS component seems to add small scale 
random noise on top of the original Gaussian CMB signal. 
We turn therefore to the use of two-point peak statistics for 
distinguising between the two maps. 

To this end, we measure the TPCF of peaks in our 
Gaussian maps, our G+S maps, and our Gaussian-GS maps. 
For each value of G/i, map size, resolution scale and beam 
size, we have generated ensembles of ~ 100 maps. Lower 
panel of Figure 2] and Figure [5] show results from averag- 
ing over 100 realizations of maps with G = 10° map at 
R = \' and $ > loo. There are obvious differences between 
the TPCF in the G and G + S maps, with the latter hav- 
ing substantially more signal on small scales. Although the 
beam erases some of this (Figures |4] and [5}, a residual effect 
remains. This signal is rather different from that measured 
in a Gaussian field which has the same Gg (what we called 
Gaussian-GS previously). So we conclude that this is indeed 
a promising method for identifying the CS component in the 
maps. 

The lower panel of Figure U shows explicitly that, al- 
though the peak counts were unable to distinguish between 
GSB and Gaussian-GSB maps (Figure 0, the TPCF on 
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Figure 4. Upper panel: Comparison of a Gaussian+String map 
(left) with a pure Gaussian map which has the same total power 
spectrum (right); blue dots show peaks above # = 0.5<7o- In both 
cases, Gfi = 8 x 10 -7 , the resolution R = V and the map size 
shown is 5° x 5°. The morphology of these two maps is quite 
different. Lower panels: TPCF of $ = I.Octo peaks in these maps 
and differences between them. 



Figure 5. Two-point correlation function of peaks above $ = 1<jo 
for different values of Gfi. Here FWHM is 4'. A£(0) corresponds 
to difference between TPCF of various cases indicated in plots. 



scales ^ 12' can. Figure [5] shows that the ability to dis- 
criminate depends on G/jl and the beam size FWHM. 



2.5 Quantitative limits 

To quantify this we first compute the Student's t-test based 



t{9) 



(5) 



where £(0) is the TPCF and cr(0) is the mean standard de- 
viation of each term in the numerator. The symbols o and 
correspond to the G + S and G measurements and to 
(G + S)B and (G)B with beam effect. For each 0, the corre- 



sponding P- value, p(0), are calculated. Degrees of freedom 
based on the t-distribution function are 2N S im — 2, where 
N S im is the number of simulated maps. 

We then define % 2 = -2^1np(0). The final P-value 
related to x 2 is calculated based on chi-square distribution 
function with 2(0 m ax — 0min)/A0 — 2 degrees of freedom. 
Fig. [6] shows this P-value as a function of Gfi for various 
maps with G = 10°. We have drawn lines at p = 0.0027, 
and 0.0455, since these correspond to 3a and 2a significance 
levels. This shows that the TPCF can detect CS at 95%CL 
provided G/jl > 1.2 x 10 -8 in maps without instrumental 
noise. If noise is present, with rmS a-noise — 10/iK, then this 
limit increases to G/jl > 9.0 x 10 -8 . Including beam smear- 
ing further degrades our limits: the minimum detectable CS 
becomes Gji > 1.6 x 10 -7 at 2a confidence interval. Table [T] 
summarizes our results. 
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Figure 6. P-value as a function of G/i. Upper left panel indicates 
the result for Gaussian CMB map. The effect of beam on the 
capability of TPCF to detect CS has been shown in upper right 
panel. Lower left panel corresponds to P-value for map in the 
presence of instrumental noise. The effect of finite beam size and 
instrumental noise have been indicated in lower right panel. To 
determine each point in this plot we did average on at least 100 
ensembles. 



3 CONCLUSION 

If they exist, cosmic strings are expected to leave an imprint 
in the CMB. We argued that although such strings may alter 
the power spectrum (Figure 1) and the statistics of hot and 
cold spots (Figure 2), the change to the one-point stastistics 
of peaks in CMB maps cannot be distinguished from that for 
a Gaussian field having the same power spectrum (Figure 3). 
On the other hand, the two-point statistics show differences 
(Figures [4] and [5]) which we believe can be used to reject the 
hypothesis that the map is a purely Gaussian (Figure [6]). 

We argued that CS will be detected at high significance 
only if the string tension is sufficiently high: G\i > 1.2 x 10 -8 . 
Accounting for the fact that instrumental noise complicates 
the measurement increases this limit to G\i > 9.0 x 10 -8 
(Figure [6] and Table [J). The CS signal is particularly strong 
on arcminute and smaller scales. Some of this signal is re- 
moved if the beam size of the experiment is larger than this 
scale. For a 4' beam, the limit is G\i > 1.2 X 10 -7 at 2a 
confidence interval. Broader beams increase the limiting G\i 
further. 

We have argued that two-point statistics of peaks (the 
pair correlation function) are better than one-point statistics 
(peak number counts) for distinguishing between models. 
Our results suggest that the n-point correlation functions of 
peaks can be used for similar purpose. This is interesting in 
view of previous work showing that the 3-point statistics of 
all pixels is not very informative. 



Map 


2a 


3(7 


G-GS 


Gii > 1.2 x 10- 8 


G/jl > 1.6 x 10~ 8 


GN-GSN 


Gfj, > 9.0 x 10- 8 


Gfi > 1.2 x 10~ 7 


GB-GSB 


Gfi > 1.2 x 10~ 7 


Gfi > 1.5 x 10~ 7 


GBN-GSBN 


Gfj, > 1.6 x 10~ 7 


Gfi > 2.2 x 10- 7 



Table 1. Limits on G\i which come from analysis of the dif- 
ferences between the TPCF for maps with and without cosmic 
strings, when R = 1', FWHM=4 / and a noise = 10 fiK. 



Final remark is that it could be interesting to use more 
realistic models ( Landriau & S hellardl 120031 ; iFraisse et al.l 
120081 : lLandriau S hellard 2010h to simulate a map taking 
all contributions of cosmic strings into account and apply our 
method to examine the effect of cosmic strings in our future 
works. Also, it is useful to apply this method to discriminate 
the effect of generic cosmic strings and superstrings. 
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